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The central black hole mass of the high-a but low-bulge-luminosity 
lenticular galaxy NGC 1332* 
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ABSTRACT 

The masses of the most massive supermassive black holes (SMBHs) predicted by the Mbh-c 
and Mbh-L relations appear to be in conflict. Which of the two relations is the more funda- 
mental one remains an open question. NGC 1332 is an excellent example that represents the 
regime of conflict. It is a massive lenticular galaxy which has a bulge with a high velocity 
dispersion a of ~ 320 km s _1 ; bulge-disc decomposition suggests that only 44% of the total 
light comes from the bulge. The M-qu-ct and the Mbh-L predictions for the central black 
hole mass of NGC 1332 differ by almost an order of magnitude. We present a stellar dynami- 
cal measurement of the SMBH mass using an axisymmetric orbit superposition method. Our 
SINFONI integral-field unit (IFU) observations of NGC 1332 resolve the SMBH's sphere of 
influence which has a diameter of ^0.76 arcsec. The a inside 0.2 arcsec reaches ~ 400 km 
s _1 . The IFU data allow us to increase the statistical significance of our results by modelling 
each of the four quadrants separately. We measure a SMBH mass of (1.45 ± 0.20) x 10 9 M Q 
with a bulge mass-to-light ratio of 7.08 ± 0.39 in the i?-band. With this mass, the SMBH of 
NGC 1332 is offset from the M-qb-L relation by a full order of magnitude but is consistent 
with the A/bh-ct relation. 

Key words: galaxies: individual: NGC 1332 - galaxies: kinematics and dynamics. 



1 INTRODUCTION 

It is now widely accepted that massive elliptical galaxies and (clas- 
sical) bulges of spiral galaxies harbour supermassive black holes 
(SMBHs) at their centres with masses ranging between ~ 10 6 - 
1O 1O M . An increasing number of SMBH detections have led to 
the discovery of empirical correlations between the SMBH mass 
Mbh and t he velocity dispersion a or the luminosity L of the 
host bulge dKormendv & Richstonel fl995 ; iMagorrian et alj|l998t 
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2002 1 ; iMarconi & Hunt! 120031 ; lHaring & Rixll2004 iGiiltekin e~ 
2009). These correlations indicate a strong connection between the 
SMBH and the host galaxy. Consequently, SMBH studies are es- 
sential for a better understanding of galaxy formation processes 
dSilk&Reeslll998h . 

Because of its crucial role, SMBH demographics has become 
one of the key ingredients in cosmological simulations and theoret- 
ical models in recent years dGranato et alj|2004l ; ISomerville et al.l 
l2008t iDi Matteo et aljHoolh . An inventory of the SMBH popula- 
tion is therefore necessary. The number of secure measurements is 
presently only ~ 50 and quiescent or weakly active SMBH beyond 
the local universe can typically not be well-studied. In this situa- 
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tion, the Mbb-o and Mbb~L relations have become valuable tools 
in predicting the SMBH mass and mass function. By combining 
the distribution of the readily observed parameter a (or L) with the 
corresponding Mbb-o (or Mbh-L) relation, the SMBH abundance 
can be indirectly estimated. 

This method, however, suffers from a bias: the Mbh-c rela- 
tion predicts fewer SMBHs with masses > 10 9 Mq than does the 
Mbh-L relation. This happens because the luminosity and velocity 
dispersion functions obtained from SDSS are differ ent from those 
in the SMBH sample which define the relati ons dBernardi et al.l 
l2007l ; lTundo et alJl2007l) . lBernardi et al.l d2007h argue that the bias 
lies in the SMBH sample. Assuming that it is due to a selection 
effect, the ir models suggest t hat the Mbh-o" relation is more fun- 
damental. lLauer et al] d2007h who examine a sample of brightest 
cluster galaxies, however, find that the local SMBH mass func- 
tion for A4bh > 3 x 1O 9 M0 inferred from the Mbh-L relation 
is in better agreement with the volume density of the most lu- 
min ous quasars. Further more, from the quasar luminosity func- 
tion, [Shjeldretal] d2006l) estimate the density of relic SMBH with 
masses > 5 x 10 9 Af Q to be 10 2 3 Gpc~ 3 . This should translate 
to the same density for galaxies with a> 500 km s" 1 based on 
the local AfsH-cr relation. No objects with such high a have been 
found in the local universe. lBernardietai] f2006) find at most two 
or three candidates in SDSS with a > 500 km s _1 in a volume of 
~ 0.5 Gpc -3 , but these might be a result of superposition effects. 
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Several other author s jNetzej 120031 ; IWvithe & Loebl 120031 ; 
IShields & Gebhardli r2Q04 ) have also noted that the a implied (via 
the Mbh-it relation) by the largest SMBH masses inferred from 
quasars exceed the largest a found in local galaxies. If we trust that 
those SMBH masses are c orrect, then the A/bh-ct relation must be 
different at the upper end. IWyithd J2006h argues that the A/bh-ct 
relation is curved upwards at the high-cr end rather than linear in 
log-log space. If this is true, then the abundance of local SMBHs 
would be closer to the quasar prediction. 

Since the shape of the A/bh-ct relation at the upper end criti- 
cally determines the space density of the most massive SMBHs, it 
is important to characterise its slope and intrinsic scatter. Currently, 
this high-cr regime is scarcely sample d and there are uncer tainties 
as to how the relation should behave. iBeifiori et all {2009) derive 
upper limits of 105 SMBH masses based on HST spectroscopy of 
ionised gas and find that the A/bh-ct re lation flattens a t the high-cr 
end, opposite to what was suggested by Wyithe (2006). To resolve 
these uncertainties, an increased sample of direct SMBH detec- 
tions in high-cr galaxies is required. If the A/bh-ct relation indeed 
breaks down in the high-cr regime, then a sufficiently large num- 
ber of SMBH mass measurements at cr~ 300 — 400 km s _1 will 
likely be able to detect this. For this reason, we undertake an ob- 
servational campaign to measure SMBH masses of high-cr galaxies 
using the near-infrared integral-field spectrograph SINFONI at the 
Very Large Telescope (VLT). 

NGC 1332 is a nearby massive SO galaxy which resides in 
the Eridanus Cloud. Its orientation is close to edge-on and the 
galaxy appears to be a normal lenticular. Adopting the bulge ve- 
locity dispersion from HyperLed^] (~ 320 km s" 1 ), we expect 
that the sphere of influence has a diameter of ~ 0.76 arcsec, 
which is resolved by our observations. The A"-band magnitude of 
NGC 1332 from the 2MASS Large Galaxy Atlas is 7.052, which 
gives a K-band luminosity of 1.56 x 10 11 Lq after a correction for 
Milky Way extinction of 0.012 magnitudes (from the NASA/IPAC 
Extragalactic Database - NED). Our photometric bulge-disc de- 
composition (Section [4. 1 b implies a bulge-to-total luminosity ra- 
tio of 0.44. Given this bulge luminosity, we would expect to 
find a M B h of 1.37 x 1O 8 M from the M B h-L k relation of 
Marconi & Hunt (2003). On the other han d, the A/bh-ct relation 
of iTremaine et al.l d2002l) or iGultekin etail d2009h constrains the 
M BH to ~ 9.0 x 10 s M e or ~ 9.7 x 10 s M respectively. The 
SMBH masses given by these two relations clearly differ by almost 
an order of magnitude which makes NGC 1332 a particularly at- 
tractive case. As a side note, the software and sources that we used 
to derive the bulge lum inosity were different from those used by 
iMarconi & Huntl d2003h . which could introduce systematic errors. 
It is, however, unlikely that the situation for NGC 1332 would be 
significantly affected. 

To date, the only published A/bh m easurement for NGC 13 32 
is provided by a recent X-ray study of iHumphrev et alj d2009h - 
hereafter H09. They make use of the Chandra X-ray Observatory 
data and rely on the assumption of hydrostatic equilibrium in the 
analysis, which results in a SMBH mass of 5.2±|;| x lO 8 A/ (with 
the Mbh-ct relation as a Bayesian prior). This estimate lies in be- 
tween the prediction of both relations, although considering the er- 
ror bars, the A/bh-ct relation is slightly favoured. It is, however, not 
clear if their black hole mass would lean more towards the A/bh- 
L relation if they were to use the Mbh-L instead of the A/bh-ct 
relation for their Bayesian prior. 
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In this paper, we measure the SMBH mass in NGC 1332 us- 
ing a stellar dynamics approach which is not biased by any of the 
two relations. Thr oughout, we adopt a distance of 22.3 Mpc from 
lTonrvetal.U200ll) . corrected for the Cepheid zer o point by apply - 
ing a distance modulus shift of -0.06 magnitudes dMei et aU2005h . 
At this distance, 1 arcsec corresponds to 0.1 1 kpc. 

This paper proceeds as follows. We describe the data and data 
reduction in Section [2] Details on the derivation of the kinematics 
follow in Section[3] The photometry, bulge-disc decomposition and 
deprojection to model the luminosity profile are described in Sec- 
tion|4] In Sections[5]and[6] we present the dynamical modelling and 
the results. Lastly, we summarise and discuss our results in Section 

m 



2 DATA 

2.1 SINFONI observations and data reduction 

The integral-field data presented here were obtained on Novem- 
ber 25, 2008, as part of the guaranteed time observations with 
SINFONI on the UT4 of the VLT. SINFONI (Spectrograph for 
INtegral Field Observations in the Near Infrare d) is a Spectrom- 
eter f or Infrared Faint Field Imaging (SPIFFI) (Eisenhauer et al. 
120031) combined with the Mu lti-Application Cur vature Adap- 
tive Optics (MACAO) module dBonnet et all |2004) . NGC 1332 
was observed in the /f-band (1.95-2.45 /im) in two different 
spatial resolutions, i.e. 0.05 x 0.1 arcsec 2 spaxel -1 (hereafter 
"lOOmas") resulting in a 3 x 3 arcsec 2 field-of-view (FOV) and 
0.125 x 0.25 arcsec 2 spaxel -1 with 8x8 arcsec 2 FOV (here- 
after "250mas"). For the former, adaptive optics (AO) correction 
was applied using the nucleus of the galaxy as the natural guide 
star. The observations followed a sequence of 10-minute exposures 
of object-sky-object-object-sky-object (O-S-O-O-S-O). Each expo- 
sure was dithered by a few spaxels to allow for bad pixel correction 
and cosmic ray removal. A total of 40 minutes on-source expo- 
sure time was obtained for the 250mas scale and 80 minutes for 
the lOOmas scale. To estimate the AO performance and the point 
spread function (PSF) due to atmospheric turbulence, we regularly 
observed a PSF star after the science exposure sequence. 

The reduction of SINFONI data was perfor med using custom 
reduc tion packages for SINFONI, i.e. ESOREX ( Modiglian i et all 
l2007t) and SPRED dSchreiber et alJl2004lAbuter et alfeOOot) . Ex- 
cept for the 250mas science data, all other frames including the 
PSF and telluric standard stars were reduced using ESOREX 
which produced cleaner spectra. Both software packages included 
all common and necessary steps to reduce three-dimensional data 
and to reconstruct a datacube. The closest sky frame was first sub- 
tracted from the science frame. The resulting frame was then flat- 
fielded, corrected for bad pixels and detector distortion and wave- 
length calibrated before the datacube was reconstructed. For tel- 
luric correction, we used three early-type stars with the spectral 
classes B3V and B5V, i.e. Hip014898, Hip023060 and Hip018926. 
As the end process, the individual science datacubes were aver- 
aged into one final three-dimensional datacube per plate scale. For a 
m ore detailed descript ion of the data reduction, we refer the reader 
to lNowak et ail d2008l) . Fig. □ shows SINFONI images of the two 
platescales which resulted from collapsing the datacubes along the 
wavelength direction. 

After the reduction, the individual PSF star images for the 
lOOmas scale were averaged, normalised and then fitted with the 
commonly-adopted double Gaussian function. To account for the 
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Figure 1. SINFONI images of NGC 1332 in two resolution scales: lOOmas 
(upper panel) and 250mas (lower panel). 



asymmetry of the PSF, both Gaussian components, i.e. a broad and 
a narrow one, were set to be non-circular. The fit gave the full width 
at half-maximum in the x and y direction (FWHM X , FWHM y ) of 
(1.03 arcsec, 0.94 arcsec) for the former and (0.15 arcsec, 0.13 arc- 
sec) for the latter. The narrow component contributed to 36% of 
the total flux and resolved the expected sphere of influence of the 
SMBH. We used this SINFONI PSF parameterisation for the sur- 
face brightness deprojection (see Section |4~3l . Fig. [2] presents the 
PSF image with the fit along the two spatial axes. We note that this 
PSF image does not strictly represent the true PSF since the ac- 
quisition was not done simultaneously with the galaxy observation. 
The fit is also admittedly not perfect, but for our purpose here the 
deviation does not lead to a significant error. We verify this and 
discuss the effect of the PSF uncertainties further in Appendix lAl 
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Figure 2. A double non-circular Gaussian fit to the SINFONI lOOmas PSF. 
The fit along the x and the y-axis of the detector are shown in the left and the 
right panel respectively. The black dashed lines are the individual Gaussians 
and the red solid lines are the sum of the Gaussian components; diamonds 
are the actual PSF. 



(first reducing the individual amplifier sections, then scaling and 
joining them into single-chip images). Since the imager was actu- 
ally a two-CCD mosaic, the result was three pairs of single-chip 
images; we combined these into a single mosaic using the SWarp 
package dBertin et al.ll2002l) . The final image had a seeing of 0.80 
arcsec FWHM (mean of Moffat profiles fitted to nine bright, unsat- 
urated stars) and a plate scale of 0.332 arcsec pixel -1 , and is shown 
in Fig.[3] 

From the HST archive we retrieved two 160s Wide-Field Plan- 
etary Camera 2 (WFPC2) images in the F8 14W filter (Proposal 
ID 599 9, PI Phillips), which were also used in lKundu & Whitmord 
d200lh . Since the images were well aligned, we combined them di- 
rectly using the STSDAS crrej task in IRAF. Finally, we generated 
A'-band images from our 250mas and lOOmas SINFONI datacubes 
by collapsing (averaging) the cubes along the spectral direction. 

We first calibrated the combined NTT-EMMI mosaic image, 
which was large enough to ensure accurate sky subtraction, using 
Cous ins R photometry from the literature ePrugniel & Heraudeaul 
1998). The WFPC2 image was then calibrated by matching surface 
brightness profiles: fitting ellipses of fixed orientation (position an- 
gle and ellipticity) to both images, then simultaneously fitting for 
the best combination of sky background (in the WFPC2 image) and 
scaling so that it matched the NTT-EMMI profile outside the central 
2 arcsec (where the worse seeing of the NTT-EMMI image affected 
the profile). Finally, we repeated the process for profiles from the 
two SINFONI images by matching them against the (scaled, sky- 
subtracted) WFPC2 profile. 



2.2 Imaging data 

We made use of three types of imaging data for NGC 1332. For 
the large-scale analysis, we searched the major telescope archives 
and found several i?-band images obtained with the red channel 
of ESO Multi-Mode Instrument (EMMI) on the 3.5m New Tech- 
nology Telescope (NTT) at La Silla. These images were originally 
taken on 2005 December 3, as part of a spectroscopic program 
(Program ID 076.B-0 182(A), PI Aragon-Salamanca). We selected 
the three best 10s exposures (the fourth exposure had strong back- 
ground variations) and reduced them with standard IRAF tasks 



3 KINEMATICS 

3.1 SINFONI kinematics 

We extracted the line-of-sight velocity distributions (LOSVDs) 
non-para metrically using a Max imum Penalised Likelihood (MPL) 
method dGebhardt et alJuOOObl) . The galaxy spectra were decon- 
volved using the weighted linear combination of a set of stellar 
templates consisting of K and M stars. These stars were previously 
observed using SINFONI with the same instrumental setups as in 
the galaxy observations. We briefly describe here the kinematics 
analysis that we have performed. It is largely similar to that in 
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Figure 3. ij-band isophotes of NGC 1332, from archival NTT-EMMI images. The left-hand panel shows the whole galaxy, plus its neighbour NGC 1331 (near 
the left-hand edge); the image has been smoothed with a 5-pixel-wide median filter. Contour levels run from 24.5 to 15.0, in steps of 0.5 mag/arcsec 2 . The 
right-hand panel shows the inner regions of the galaxy, including the bulge, using an unsmoofhed image; contour levels run from 23.1 to 15.0, in steps of 0.3 
mag/arcsec 2 . North is up and east is to the left. 



iNowak et al.1 J2007I) and lNowak etal] d2008h . so we refer the reader 
to those papers for details on the kinematics derivation. 

Reliable kinematics from MPL can be obtained when the 
S/N is sufficiently high. Therefore, to optimise and homogenise 
the S/N, we binned th e pixels into angular and radial bins as in 
iGebhardt et al] J2003I) by luminosity-weighted averaging of the 
spectra. The galaxy is divided into four quadrants bordered by the 
major and minor axes. Each quadrant is divided further into five 
angular bins. The centres of those bins are at the angles of 5.8°, 
17.6°, 30.2°, 45.0° and 71.6°. For the lOOmas data, there are seven 
or eight radial bins while for the 250mas data, 12-13 of those bins 
were needed to cover the FOV. We then performed the MPL method 
on the binned spectra as follows. 

We first normalised the galaxy and stellar template spectra by 
dividing out the continua. The combined stellar template was con- 
volved with a binned initial LOSVD. The LOSVD and the weights 
of the templates were iteratively changed until the convolved com- 
bined spectrum matched the galaxy spectrum. This fit was done by 
minimising the penalised \ 2 '- Xp = X 2 + a P- A certain level of 
smoothing was applied to the LOSVD via the second term where 
P, the penalty function, is the integral of squared second deriva- 
tive of the LOSVD. The smoothing parameter a determines the 
level of regularisation and its value depends on the velocity dis- 
persion of the galaxy and the S/N of the data. We estimated the 
appropriate smoothing for our data from the kinematic analysis of 
a large dataset of model galaxy spectra. These models were created 
by broadening the template spectrum with a velocity dispersion of 
400 km s~ 1 . Our data have a high S/N which reaches 70 in the cen- 
tral pixel. After the binning, the S/N increased to ~ 90 (lOOmas) 
and ~ 83 (250mas) on average. With those S/N values, the appro- 
priate values for a found in the above simulations are on average 
~ 5 and ~ 6, respectively. 

To derive the LOSVDs, we specifically fitted the first two CO 
bandheads CO(2-0) and CO(3-l) in our spectra. To minimise the 
error due to template mismatch, we measured the equivale nt width 
of the first CO bandhead as in ISilge & Gebhardl 020031) and se- 
lected only stars with similar equivalent widths for the templates. 



Across the FOV of SINFONI 250mas, the measured values range 
from ~ 11 to ~ 15A. 

We calculated uncertainties for each LOSVD from 100 
Monte Carlo realisations of the galaxy spectra. These spectra 
were obtained by convolving the measured LOSVD with the 
stellar templates. Each spectrum differs from the others in the 
amount of Gaussian noise added to the spectrum at each wave- 
length position. From every spectrum, an LOSVD was extracted 
and used to estimate the errors. For illustration purposes, we 
parametrised the LOSVDs in terms of Ga uss-Hermite moments 
(van derMarel&Franxlll993l : lGerharj[l993l) . i.e. velocity v, ve- 
locity dispersion a and two higher order terms which measure the 
asymmetric and symmetric departure from a pure gaussian veloc- 
ity profile /13 and Ha. The typical errors that were derived for v and 
a are 8.45 km s _1 and 8.87 km s -1 (lOOmas), respectively; for 
250mas, the errors are 7.92 km s _1 and 9.21 km s _1 . For /13 and 
I14, the errors are typically 0.02 for both scales. 

We present kinematic maps of NGC 1332 for both scales in 
Fig- El A significant rotation is shown by the well-ordered pattern 
in the velocity map and the anti-correlating I13. The velocity dis- 
persion is peaked at around 400 km s~ and there is a rather steep 
decline towards the outskirts. In spite of the presence of dust in the 
nucleus, the kinematic centre seems to coincide with the photomet- 
ric centre. 



3.2 Long-slit kinematics 

Our SINFONI data provide the required resolution to allow for 
an accurate measurement of the SMBH mass at the expense of 
the FOV size. To obtain constraints on the orbital structure at 
larger radii, we uti lised long-slit data which were reported in 
iKmiken et all i 19961) - hereafter K96. They observed NGC 1332 
using the Red Channel Spectrograph at the Multiple Mirror Tele- 
scope. The optical data were taken along the major axis of NGC 
1332 through a 1.25 arcsecx3 arcmin slit with a spectral resolu- 
tion of 2.6A (instrumental <r=63 km s _1 ). The major axis position 
angle (PA) given in K96 is 148° following RC3. This is about 30° 
higher than the PA we measured from the WFPC2 image, i.e. 117°. 
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Figure 4. Kinematics of NGC 1332 derived from SINFONI data, illustrated in terms of Gauss-Hermite moments (v, a, /13 and /14) as written on the side of 
each map. The first row displays lOOmas data while the second row shows 250mas scale kinematics in which the spatial scope of the lOOmas data is outlined. 
The map shows the division into radial and angular bins as described in the text. The major and minor axes of the galaxy are aligned with the abscissa and the 
ordinate of the coordinate system respectively. In the first map, the numbers (1-4) correspond to the quadrant numbers and the arrows indicate the orientation 
to the north (long) and the east (short). 



The latter PA is consistent with our SINFONI data, and also with 
the PA given by HyperLeda and ESO/Uppsala Survey. By simply 
looking at images of NGC 1332, a PA of 148° can be easily re- 
jected. We have confirmed that the slit was indeed placed along the 
long axis of the galaxy and therefore the quoted PA in K96 is a typo 
(K. Kuijken, private communication). 

The velocity profile of the slit data was derived usin g 
the Fourier Correlation Quotient (FCQ) method teendej [l990h . 
parametrised into Gauss-Hermite parameters v, a, hz, /14. The 
LOSVDs were then reconstructed from those moments as it is the 
full LOSVD that is fitted in the modelling. The data provide kine- 
matics up to a radius of ~ 80 arcsec on both sides of the major 
axis. 

The SINFONI and the long-slit datasets overlap within 3.5 
arcsec and they are broadly consistent with each other at radii > 2 
arcsec. Inside that radius, v and a derived from the long-slit data are 
systematically lower. For instance, at ~ 0.2 arcsec, the slit data give 
a a of ~ 335 km s _1 while the 250mas and the lOOmas SINFONI 
data give ~ 365 km s^ 1 and ~ 385 km s _1 respectively. Moreover, 
the SINFONI a also shows a steeper gradient. We suspect that these 
differences are due to the seeing during the slit observation and the 
better resolution of the SINFONI data. Since the PSF of the slit 
observation is not known, it is difficult to investigate the discrep- 
ancy. However, we find that within the discrepant region SINFONI 
kinematics dominates and renders the slit kinematics unimportant, 
as briefly explained in the following. We modelled NGC 1332, as 
described in Section[5] using both datasets simultaneously. We ran 
several models using identical set-ups, varying only the slit PSF 
to several reasonable values (1.5 arcsec and 2 arcsec). We also ran 
models where we used only slit data outside the SINFONI FOV 
(> 3.5 arcsec). We are able to show that all those variations do 
not alter the modelling results. With this finding, whenever we in- 
clude slit kinematics in the dynamical modelling, we consider only 
slit datapoints outside the SINFONI FOV to ease the computational 
load and time. 



4 STRUCTURAL AND LUMINOSITY MODELLING 



4.1 Bulge-disc decomposition 

Since NGC 1332 is an SO galaxy with a fairly well-defined bulge 
and disc, we investigated bulge-disc decompositions, both for mod- 
elling purposes (i.e., in case the bulge and disc stellar populations 
might have different mass-to-light ratios) and so that we could esti- 
mate its bulge luminosity in order to see where NGC 1332's SMBH 
fell in the SMBH-bulge relations. We tried two approaches, which 
agreed extremely well. The first involved fitting free ellipses (i.e., 
with position angle and ellipticity allowed to vary) to the isophotes 
of the NTT-EMMI and SINFONI images (carefully masking out 
the dust lane in the SINFONI images) and then combining these 
into a single 1-D surface-brightness profile (see Section l4~2"t . The 
best fit to this profile was a combination of an exponential (repre- 
senting the outer disc, with central surface brightness /Lto = 18.78 
and scale length h — 33.6 arcsec), a Sersic function (represent- 
ing the bulge, with n = 2.36, surface brightness /i e = 18.36, and 
effective radius r e = 9.15 arcsec), and a small, narrow Gaussian 
representing a possible central star cluster. 

We also tried a 2D bulge-disc deco mposition, using ver- 
sion 2.2 of th e BUDDA image-fitting code jde Souza et al.ll2004l : 
Gadotti 2008) and the NTT-EMMI image by itself; we masked the 
dust-affected part of the nucleus and turned on seeing correction. 
The best fit was with an exponential disc having an ellipticity of 
0.73 ([io = 18.66, h — 32.96 arcsec), a Sersic bulge with elliptic- 
ity = 0.27 (n = 2.34, /i e = 18.18, r e = 8.39 arcsec), and a small 
point source (representing, e.g., a stellar nucleus) contributing to 
1.2% of the total light. The agreement with the 1-D decomposition 
is excellent. We note that the 2D decomposition corresponds to a 
bulge-to-total light ratio of 0.43 (0.44 if we include the small point 
source), so NGC 1332 is, despite its high central velocity disper- 
sion, still a (marginally) disc-dominated galaxy. 
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4.2 Photometric models 

We constructed two main photometric models. The first was a 
single-component model, where the galaxy was represented by 
a single surface-brightness profile, with variable ellipticity an d 
higher-order moments £14, a^, etc. dBender & Moellenhofd[l987h . 
To construct this profile we fit ellipses to isophotes using three im- 
ages: the NTT-EMMI J?-band image and the two SINFONI im- 
ages (i.e., lOOmas and 250mas datacubes collapsed along the wave- 
length axis). For both SINFONI images, we masked out the regions 
affected by the nuclear dust lane before fitting ellipses. The NTT- 
EMMI profile was used for r > 4.5 arcsec; at smaller radii, where 
better resolution and lower dust extinction were needed, we used 
the SINFONI profiles (the 250mas profile was used from 4.5 arc- 
sec in to 1.5 arcsec, with the profile at r < 1.5 arcsec coming from 
the lOOmas image). We see no noticeable gradient in V — I colour 
from HST images in the innermost 13 arcsec along the south semi- 
minor axis, which is less affected by dust. Since dust regions were 
masked out during profile extraction, V — I should be flat, and so 
we do not expect a R — K color variation, either. The combination 
of images from the R and K bands can therefore be justified. 

The second photometric model involved separating the galaxy 
into a bulge and a disc, with potentially different stellar mass-to- 
light ratios. Since the 2D bulge-disc decomposition worked well, 
we decided to use the disc component from that fit as the disc com- 
ponent for modeling purposes. Because we wanted to match the ac- 
tual surface brightness and isophote shapes in the central regions as 
closely as possible, using the (fixed-ellipticity) Sersic bulge model 
from the 2D fit would have been too crude an approximation; even 
including the small Gaussian component still produces residuals 
and does not reproduce the actual ellipticity profile of the galaxy's 
inner regions. Instead, we first subtracted the disc model from our 
images (specifically, the NTT-EMMI image and the two SINFONI 
images) to create residual, "bulge-only" images and then performed 
variable-ellipticity fits on these images in order to generate a bulge 
model for modeling purposes. The result of this was that the surface 
brightness and ellipticity of the bulge component tracked the actual 
brightness and isophote shapes in the very central regions as ac- 
curately as possible; moreover, the combined light from this bulge 
component and the 2D disc model reproduces the original galaxy 
light distribution. 

To construct the final bulge component, the ellipse fits to the 
NTT-EMMI image were used for r > 4.5 arcsec; the 250mas SIN- 
FONI image was used for r = 1-4.5 arcsec, and the lOOmas im- 
age was used for r < 1 arcsec. As was the case for the single- 
component model, we carefully masked out the dust lane in both 
SINFONI images before running the ellipse-fitting software. At 
radii > 15 arcsec, the NTT-EMMI profile from the residual image 
became significantly affected by deviations of the disc from the 2D 
model which we subtracted from the image. To ensure a relatively 
smooth luminosity model, we replaced the data at r > 15 arcsec 
with an extrapolation of the best-fitting Sersic component from the 
2D fits ( section |4~Tt , fixing its ellipticity to 0.27. At these radii, the 
light is dominated by the disc component, so small variations in the 
bulge component have minimal effect on the modelling. 

Finally, we constructed an alternate single-component profile 
for testing purposes (see Appendix IAt : this was identical to the 
SINFONI+NTT-EMMI profile described above except that ellipse 
fits to the WFPC2 F814W image were used for r < 0.5 arcsec. 
Problems with strong dust extinction in the WFPC2 image ulti- 
mately led us to reject using this data for the actual modelling. 

Isophotal shapes of all the models (the one-component model, 
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Figure 5. The isophotal shape analysis of NGC 1332. As a function of 
the logarithm of the semimajor axis distance we show the .R-band surface 
brightness (top left), the ellipticity (top right), the <I2 (middle left), 04 (mid- 
dle right), ae (bottom left) and ag (bottom right) coefficients of the isopho- 
tal Fourier analysis. The solid lines refer to the one-component photometry. 
The short dashed lines show the photometric parameters as derived from the 
HST image in the inner 0.5 arcsec. The dotted and long dashed lines show 
the bulge and disc components of the two-component photometry, respec- 
tively. 



bulge and disc of the two-component model and the single- 
component model using WFPC2/HST image) is shown in Fig. [5] 



4.3 Deprojection 

The dynamical modelling requires knowledge of the three- 
dimensional luminosity distribution v (cf. Section[5j. Each photo- 
metric profile, i.e., the single-component model, the bulge and the 
disc model was deprojected separately, resulting in the (i?-band) 
3D luminosity profiles v, V\, and Vd respectively. Since the disc 
flattening is small, the orientation of NGC 1332 must be close to 
edge-on. We therefore assume an inclination of 90° throughout the 
paper unless otherwise stated. The edge-on deprojection gives an 
intrinsic flattening of 0.27 for the disc. 

The deprojections of the bulge, the disc and the one- 
component model were done under the ass umption of axisy mme- 
try using a modified version of the code o f Mag orriarJ <1999h . We 
briefly describe it as follows. We firs t constructed an initial den- 
sity model, defined in equation 10 of iMagorriarJ (1999). This was 
done by running through a grid of the model parameters. Each of 
the models was convolved and projected to be compared with the 
observed surface brightness. The model with the smallest x 2 was 
selected as the initial model. This model was then refined by ap- 
pl ying small changes through a simulated annealing procedure as 
in Magorrian (1999). The convolution, projection and comparison 
steps were subsequently repeated after each change. The iteration 
was stopped and the final density model was reached when the 
model matched the observations within a pre-determined accuracy. 
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Figure 6. The model of PSF-deconvolved luminosity density of the one- 
component (red dashed line) and the two-component (black solid line) den- 
sity models for an inclination of 90°. The subscript "R" refers to the R- 
band. Thick and thin lines refer to density profiles along the major and mi- 
nor axes respectively. 

For the disc component, the PSF effect is negligible because 
the disc is very faint at the innermost radii where the PSF becomes 
important (the central luminosity density of the disc is at least three 
orders of magnitude lower than that of the bulge). For the bulge 
and the one-component model, where the innermost isophotes are 
based on the SINFONI lOOmas image, we had to take the PSF 
into account during the deprojection. For this purpose, we used the 
double-Gaussian parameterisation of the SINFONI PSF described 
in Section [2~Tl to implement the PSF convolution. 

We show the PSF-deconvolved luminosity density profiles for 
both density models in Fig. [6] For the two-component models, we 
plot the total quantity (bulge + disc). The density of the disc ex- 
ceeds that of the bulge at r > 12 arcsec. The difference of the den- 
sity profiles at < 1 arcsec reflects a certain level of uncertainties in 
the photometry and the bulge-disc decomposition, which gets am- 
plified in the deprojection. The luminosity densities differ mainly 
along the minor axis. It is apparent that the two-component density 
model is rounder in the centre. We investigate how these uncertain- 
ties and differences affect the SMBH mass estimate in Section[6] 



5 DYNAMICAL MODELLING 

We em ployed an orbital synthesis method based on Schwa rzschildl 
dl979t) to model NGC 1332 and to measure the SMBH mass. 
In particu lar, we used the three-integr a l axisymmetric code de- 
scribed inlGebhardt et al] d2000bLl20rB) . |Thomas etail J2004) and 
Siopis et al.l d2009h . The modelling procedure includes the follow- 
ing steps. (1) Calculation of atrial gravitational potential consisting 
of the contribution from stars and the SMBH. (2) Generation of an 
orbit library which obeys the given potential. (3) Calculation of 
orbital weights such that the orbit superposition satisfies the light 
distribution and the kinematical constraints. (4) Repetition of steps 
(l)-(3) for different trial potentials obtained by varying the mass- 
to-light ratio(s) and the SMBH mass. The best-fitting parameters 
are then found through a \ 2 analysis. 

We used each of the two density models described in Sec- 
tion [43] for the modelling. For the one-component density model, 
the mass distribution of NGC 1332 followed p = Tv+M-Bii8(r), 
where the stellar mass-to-light ratio T and A/bh were the only free 
parameters. For the two-component density models, the mass dis- 
tribution became p — YbVb + YdVd+MBnS(r) with the bulge 



mass-to-light ratio Yb, the disc mass-to-light ratio Yd and Mbh as 
the free parameters. 

Each of our orbit libraries consisted of about 2 x 15000 orbits, 
i.e. two identical sets of 15000 orbits, opposite in angular momen- 
tum directions. We set the maximum radius of the libraries to 100 
arcsec and limit the extension of the long-slit data to 30 arcsec. 
For each of the modelling runs that we conducted (cf. Section [6J, 
we modelled each quadrant separately, resulting in four different 
SMBH mass values, one for each quadrant. This highlights the ben- 
efit of having integral-field data, i.e. (1) to justify the assumption of 
axisymetry, (2) to have a complete spatial coverage of kinematic 
data and (3) to allow for four independent measurements of the 
SMBH mass (under the assumption of axial symmetry). In order 
to implement the PSF convolution in the modelling as accurately 
as possible, we directly used the SINFONI PSF images from both 
scales, instead of the double Gaussian parameterisation. 

The SINFONI data were mapped into 226 spatial bins (see 
Fig. |D, 120 of which come from the lOOmas data and sample the 
inner part. As kinematic constraints, we used the full LOSVD of 
each bin, sampled into 25 velocities from -1543 km s _1 to 1543 
km s _1 . These provided a total of 226 x 25 = 5650 kinematic 
observables to be fitted by the dynamical models (about 1412 ob- 
servables per quadrant). The LOSVDs from the long slit data (3.5 
arcsec< r <30 arcsec) were binned in the same way, increasing 
the total number of observables by 650 (~ 162 for each quadrant). 



6 RESULTS 

6.1 The black hole mass and stellar mass-to-light ratio 

As described in the previous sections, we prepared two kinematics 
datasets and constructed two stellar density models. For the latter, 
an inclination of 90° is assumed. In this subsection we discuss the 
results for the edge-on models. We address the uncertainties due to 
inclination in another subsection i6.3l , where we also show that the 
edge-on models produce better fits to the data. 

We performed five different modelling runs using different 
combinations of the density models and kinematics to probe the 
possible systematic uncertainties that could arise from various con- 
straints/datasets. For each density model, there were two runs, one 
with and one without slit kinematics (lOOmas and 250mas SIN- 
FONI kinematics were always used). We list all runs with the re- 
sulting Mbh for each quadrant in Table Q] and similarly for the 
mass-to-light ratios in Table [2] The naming of the five runs is de- 
scribed in the caption of the first table. Whenever the slit kine- 
matics was used, only slit datapoints within 3.5 arcsec< r <30 
arcsec were included in the modelling. For the runs with the two- 
component density model, we first used identical values for Tf, and 
Yd (runs 2A and 2B) and then we let both parameters vary with 
respect to each other (run 2B*). For run 2A, there were 14 values 
of M BH (5 x 10 s M Q to 3.6 x 10 9 M Q ) and 20 Y b (3 to 13). In 
run 2B, we calculated models for the same 14 values of Mbh as in 
run 2A and 12 Y b (4 to 10). A set of 14 T d (2 to 1 1) values were 
then added for run 2B*. Run 1 A and IB used 20 Mbh (5 x 10 8 to 
5 x 1O 9 M ) and 20 T (1 to 20). 

The corresponding \ 2 curve for Table 1 is shown in Fig. [7J 
We plot the normalised % 2 for all the models and all the quadrants. 
Ideally, the normalisation is done by dividing the x 2 by the number 
of degrees of freedom (dof). For the SINFONI data, the smoothing 
parameter included in the LOSVD derivation correlates the veloc- 
ity bins, such that the number of dof is smaller than the number of 
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Table 1. The best-fitting Mbh, marginalised over all mass-to-light ratios, for the four quadrants and their average for all runs that we performed. All values 
are stated in units of 1O 9 M . The lcr errors (~ 68% confidence level), derived from the A\ 2 analysis, are given for each quadrant. The last row gives the 
average of the measurements of the four quadrants with the average of their 1-cr errors. We compare these errors with the parenthesised values next to them 
which are the standard deviation (rms) of the best-fitting values of the four quadrants. We expect them to be comparable (see Section [6}. The naming of the 
runs is chosen as follows. The letters A and B identify the datasets used for the kinematics. Run A used SINFONI data only while run B also used slit data 
within 3.5" < r < 30". The numbers 1 and 2 indicate the use of the one-component and two-component density models respectively. In ran 2A and 2B we 
set Tt,=T [ j, whereas in ran 2B* we decoupled Tj, and Tj. All runs adopted an inclination angle of 90°. 



Run 1A Run IB Run 2A Run 2B Run 2B* 

(single-component, (single-component, (two-component (Tj,=T,j), (two-component (Y^Y^), (two-component (Y^Yj), 





SINFONI) 


SINFONI+slit) 


SINFONI) 


SINFONI+slit) 


SINFONI+slit) 


Quadrant 1 
Quadrant 2 
Quadrant 3 
Quadrant 4 


nq7 +0.54 
u - 9 '-0.06 
i ni+0.05 
L - ZL -0.4,1 
() 07+0-34 

1 45+0-07 


i ir+0.19 
1 ' 40 -0.13 

97+ ' 19 

1 91+O.2I 

^-o.os 

1 91+O.26 
1 - zl -0.04 


1 91+0.31 
1 - Z1 -0.12 
1 91+0.23 
^^-O.lS 
1 45+0-07 

1 45+0-12 


1.681°;°? 
1 91+0.28 

1 ' zl -0.08 
1 4^+0-17 
i - 4o -0.07 
1 45+017 


, co+0.09 
1-O«_ 26 

1 91+O.23 
J " ZA -0.08 
1 45+0-17 

1 45+021 


Mean 


1.15t°;|(0.23) 


1.21+^(0.20) 


1.33±3:if (0.14) 


1.45l{j;^(0.20) 


1.46±g;iS(0.20) 



Table 2. The best-fitting T or T D in the ij-band, marginalised over all Mbh, for the four quadrants and their average for all runs listed in Table[T] The lcr 
eiTors (~ 68% confidence level), derived from the A\ 2 analysis, are given for each quadrant. The last row gives the average of the measurements of the four 
quadrants with the average of their lcr errors. Descriptions of the runs are given in the caption of Table[T] For run 2B* where Tj, and Y^ were decoupled, we 
write down only Y since Y^ was poorly constrained. The four measurements of Y^ in ran 2B* fell within a range of 5.0 to 9.0 with an average of 8.0 and 
rms of 2.0 





Run 1A (Y) 


Run IB (Y) 


Run 2A (Y 6 =Y d ) 


Run 2B (Y 6 =Y d ) 


Run 2B* (Y b ) 


Quadrant 1 
Quadrant 2 
Quadrant 3 
Quadrant 4 


9-581818 
7.68l8;ll 
10.53i?-j! 

7 fi o+0.92 
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7.87l8;||(0.26) 


7.08l8;l(0.26) 


7.08l8;3|(0.26) 



observables (more details can be found in iGebhardt et alhoOObh . 
Because the exact number of dof is unkown, we approximate it by 
the number of observables when normalising the chisquare for runs 
1A and 2A. For runs IB, 2B and 2B*, there are an additional dof 
contributed by the slit data. Each of the slit LOS VDs was generated 
from four Gauss-Hermite parameters (see Section POt so there are 
four dof for every slit LOSVD. The normalisation factor for each 
run and quadrant is included in Fig. [7] The normalised \ 2 values 
are less than unity, as is commonly found in the orbit-based dynam- 
ical modelling. This is partly because the number of observables is 
larger than the effective number of dof due to the aforementioned 
smoothing parameter. 

For each marginalised \ 2 curve of Mbh (see Fig. [7} or Y we 
approximated the lcr error (a A\ 2 of 1) by a polynomial fitting. To 
check whether the measurements in the four different quadrants are 
in agreement with each other, we compare the lcr errors with the 
standard deviation/rms derived from the four measurements. We 
expect the rms to be comparable to the lcr errors as the latter in- 
dicate the range of values within which the measurements (from 
quadrant to quadrant) would fluctuate. For most of the results in 
Table 1 and 2, the rms do not fall far from the la errors, especially 
for A/bh measurements. It follows that (1) the four quadrants give 
consistent results and (2) the x 2 analysis provides a reliable mea- 
surement of the errors in each quadrant. 

All models strongly suggest the presence of a central black 
hole with a mass of at least 1O 9 M0. No-SMBH models were not 
included in any of the runs as for these models, A\ 2 > 20 (com- 
pared to the best-fit model) can be readily inferred from the \ 



curves. The signature of the SMBH is strongest in the models that 
include long-slit data, i.e. they give higher best-fit Mbh in most 
cases and exclude the Mbh=5 x 10 s Mq with higher confidence 
(see Fig. [TJ. M B h=5 x 10 s Mq is the lowest Mbh that we mod- 
elled and it coincides with the resulting A^bh derived from X-ray 
data by H09. It is, however, not favoured by any of our runs. The 
A\ 2 of the best-fitting model for Mbh=5 x 10 8 Mq is larger than 
10 for run 1 A and larger than 20 for the other runs. 

From the result of the five runs, we see that Mbh is rather 
sensitive to the change of the photometry and therefore th e pre- 
scribed density profiles, as is also found bv lSiopis et ail d2009l) . As 
the galaxy becomes rounder in the inner part of the two-component 
model, Mbh becomes systematically higher and Y becomes lower 
compared to the one-component model. The mass-to-light ratio dif- 
ference between runs 1A and 2A is larger than the difference be- 
tween runs IB and 2B. The addition of slit data reduces the dif- 
ference in mass-to-light ratio but the trend remains. The increase 
of Mbh in two-component models is probably related to the flatter 
slope of the density inside ~ 0.3 arcsec. It is, however, reassuring 
to see that the differences in the two mass models do not push Mbh 
beyond their lcr errors. 

The addition of the slit data seems to bring Yt to a lower 
value which ultimately increases Mbh to preserve the enclosed 
mass. Another important aspect from including the slit data is that 
it should reduce the statistical uncertainties in the modelling. This 
is indeed the case as can be seen in Table 1 and 2: lcr errors de- 
crease from run 1A to IB and from run 2A to 2B. Fig.|7jprovides a 
more straightforward way to assess this. The \ 2 curves of the run 
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Figure 7. The normalised \ 2 vs Mbh marginalised over all mass-to-light ratios for the runs listed in Table[T]for four different quadrants. The normalisation 
factor, by which we divide x 2 > is written for the individual runs in each quadrant. The best-fitting A/bh is given by the model with minimum normalised \ 2 
in the corresponding run and quadrant. These SMBH masses can be found in Table[T] 



IB are narrower and is enveloped by run 1A. The same is true for 
run 2 A with respect to run 2B and also run 2B*. 

In run 2B*, we repeated run 2B but allow for T& and Yd to be 
different from each other. The resulting A/bh and Yt in all quad- 
rants are unchanged with respect to run 2B. The best-fitting models 
of the four quadrants in run 2B* produce generally better \ 2 values 
since they were given extra freedom to fit the data. It appears that 
untying Yd and T;, does not lead to any change (or small if any- 
thing at all) in A/bh- Yd in run 2B* is not constrained well by the 
data since the kinematics data extend only up to 25-30 arcsec, ap- 
proximately where the disc becomes important. The results suggest 
that the disc has a higher mass-to-light ratio than the bulge (see the 
caption of Table[2]l. Although we varied Yd quite extensively, T;, 
is stable, showing that Yd is not correlated with T^. 

In the dynamical modelling, the resulting mass-to-light ratio is 
closely connected to A/bh • It is therefore important to constrain the 
orbital structure in the outer part as much as possible; this gives a 
stronger preference for the runs which include slit data. The result- 
ing Mbh of the three runs with slit data (IB, 2B and 2B*) are all 
consistent with each other within their la errors. As for the density 
model choice, we are more inclined towards the two-component 
density model. The bulge-disc decomposition was done based on 
the morphology of the galaxy and to allow for different components 
to have stellar populations with potentially different mass-to-light 
ratios. The results of run 2B* hint that Td is higher than Y&, al- 
though the former is not well-constrained by the data. It turns out 
that even though Y;, and Yd might be different, A/bh does not 
depend on a possible Y gradient outside the bulge. Nevertheless, 
because in run 2B* we have explored the most degrees of freedom 
we consider its Mbh as the least biased estimate compared to the 
other runs. We selected the best-fitting model in run 2B* as our 
preferred model and quote the mean of A/bh and Tt estimates in 
four quadrants as our best estimate. To be conservative, we adopt 
the largest nominal of errors as the final error margin, i.e. 0.20 for 
the Mbh and 0.39 for the Y b . 

The major-axis kinematics and the model fitting are presented 
for the four quadrants in the first and second row of Fig. [8] We 
plot v and a to represent the kinematics. We stress that it is the 
LOSVD (derived from the SINFONI spectra and from the Gauss- 
Hermite parameters of the slit data) that is fitted by the model, not 
the velocity moments. Our preferred model is shown by the red 
line. For comparison, we overplot the best-fitting model for Mbh 



of 5 x 10 s Mq (blue line). In the centre, the blue line falls below 
the red line as expected from the SMBH masses of the models. It 
climbs over the red line at intermediate radii and lies above the red 
line in the outer region because Y(, of the best-fitting model for 
Mbh of 5 x 1O 8 M0 is higher than that of our preferred model. 
This effect is shown by the dashed green line which represents the 
models where A/bh is 5 x 10 8 and T b and Y d are the same as those 
of the preferred model. At a first glance, it might not be directly 
obvious which model fits the data best. To investigate this, we plot 
the x 2 differences between the preferred model and the two models 
with A/bh =5 x 1O 8 M0 (averaged over the five angular bins) as a 
function of radius in Fig. [8] 

The A\ 2 between our preferred model and the best-fitting 
model with M B h= 5 x lO 8 A/ is the solid line while the A\ 2 
between the preferred model and the other model is the dashed line. 
For both cases Ax 2 =Xpreferred — xle8- Both lines generally lie at 
Ax 2 < 0. It is clear that models with M B h=5 x 1O 8 M produce 
worse fits than the preferred model. There seems to be a trade-off 
and an inconsistency in the models with A/bh=5 x 10 s Mq with 
respect to the data, i.e., the mass-to-light ratio that better fits the slit 
data yields a worse fit to the SINFONI data. The importance of the 
slit data in constraining the mass-to-light ratio is apparent where 
the dashed lines lie mostly above the solid lines in the region cov- 
ered by the slit data. For the solid line, A\ 2 is slightly lower in this 
region than in the SINFONI region, giving the impression that the 
slit data plays a more significant role in determining the black hole 
mass through the mass-to-light ratio. However, the dashed lines 
show that when the mass-to-light ratios were held fixed between 
models with different A/bh, there is still a strong preference for 
the preferred model from the SINFONI data. Both datasets have 
their own contributions in the SMBH mass determination and they 
complement each other. 

One might argue that the presence of dust in this galaxy would 
bias the kinematics, mass tracers and the SMBH mass estimate. We 
compared the kinematics extracted from the dust-extincted and the 
dust-free part of the galaxy and we did not find systematic differ- 
ences. There is also no asymmetric pattern in the SINFONI kine- 
matics map that can be attributed to the dust (compare Fig. Q] and 
Fig. EJl. The dust seems to be confined to a region inside ~ 0.5 
arcsec and this was masked out in the photometric analysis for 
the mass models. We do not observe any large-scale dust struc- 
ture in NGC 1332 that would significantly affect our results (cf. 
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Figure 8. Top and middle rows: kinematics fit of the models to the data along the major axis, illustrated in terms of v and <r. Circles, diamonds and squares 
represent the SINFONI lOOmas, 250mas and the long slit data, respectively, with their corresponding errors. Our preferred model from run 2B* is shown by 
the red line. The blue line shows the best-fitting (smallest x 2 ) model for Mbh=5 x 10 s Mq. The green line represents the model with A/bh=5 X 10 8 Mq 
and the same Tf,, as in our preferred model. Bottom row: x 2 differences between the three models shown in the top two rows, averaged over all angular 



bins in each quadrant. Following the colour coding above, solid and dashed lines represent Xr e( j — Xbi ue and Xred 
dash-three-dotted line is plotted along Ax 2 =0 to guide the eye and the vertical lines separate the different datasets. 



respectively. The horizontal 



iBaes et all bOOd) for a discussion of the dust effects). A minor 
influence would probably lower SMBH mass slightly, due to the 
underestimated Tf,. 

6.2 Comparison with stellar population models 

We compared the dynamical mass-to-light ratio to that of a single 
stellar population (Tssp). Estimates of stellar population proper- 
ties (age and metallici ty for the purpose of deriving Tssp) of NGC 
1332 can be found in lBarr et alj J2007h and H09. The former mea- 
surements were based on spectra from NTT-EMMI observations 
while t he latter made us e of Lick indices provided in Trag er et"al] 
Jl998h . lBarr et al.1 J2007h quote a log age of 1.103 ± 0.024 (age of 
12.7±0.7 Gyr) and a metallicity of 0.270±0.023 while in H09, the 
age is 4zti"t Gyr with a metallicity of 0.32 ± 0.3. The huge differ- 
ence in the ages (apart from the errors) result in very different val - 
ues of Tssp, as derived from the SSP models of Mar astonl d2005h . 
Compared to Tssp, our mass-to-light ratio is almost two times 
larger when using the age of H09 with a Salpeter IMF (Tssp~ 3.6) 
or about three times larger with a Kroupa IMF (Tssp~ 2.3). 

Our dynamical mass-to-light ratio is much more consistent 
with Tssp when the galaxy age is taken to be ~ 12 Gyr. This age 
is the one measured by Barr et al. and it is compatible with the up- 
per Jirmt_ofIiro rneasurement. In addition, Lick indices measured 
bv lQgandoeTaLl j2008l) also imply an old age of at least 12 Gyr. 
The Barr et al. age results in Tssp of ~ 5.0 for a Kroupa IMF or 
~ 7.9 for a Salpeter IMF. Our dynamical mass-to-light ratio falls in 
between these two values, but is closer to the Salpeter one. The ten- 
dency that the dynamical T agrees better with the Salpeter-based 
Tssp, rather than Kroupa, in mass i ve ear ly-type galaxies was pre- 
viously found by Cappella ri et alj d2006h and also by Thomas et 



al. (in preparation). Although not conclusive, this might indicate 
that the Salpeter IMF is a better represent ion of stellar pop ula- 
tions in massive early-type galaxies (see also lGrillo et al ]|2009l and 
iTreu et alj|2009h - 



6.3 Inclination effects 

According to its high ellipticity (reaches ~0.7), the orientation of 
NGC 1332 must be close to edge-on. To check for any influence 
of the residual small uncertainty in the inclination, we ran models 
with i < 90°. Constraining the disc to an intrinsic flattening no 
less than 0.2 gives 80° as the lower limit of the inclination. Using 
this angle, we repeated the deprojection step and the modelling for 
each of the four quadrants as in run 2B*. We found A/bh=L45 x 
10 9 M Q (rms=0.20 x 10 9 M Q ), T(,=6.95 (rms=0.30) and T d =7.0 
(rms=1.83). Compared to run 2B*, A/bh does not change and there 
is only a slight decrease in both T;, and T^. The uncertainties that 
could arise from the inclination assumption appear to be negligible 
in the case of NGC 1332. Compared to the preferred model (i — 
90°), the best-fitting model with i — 80° is worse by a Ax 2 ~30- 
50 in every quadrant. 



6.4 Dark matter halo 

Recently, concerns have emerged that dynamical models which do 
not include a dark halo component under estimate the true Mbh as 
is the case in M87 dGebhardt & Thomasl l2009). The reason is that 
models without dark halo require a higher stellar mass-to-light ra- 
tio in the outer part to compensate for the missing dark mass. In a 
single-component (density) model this leads to an overestimation 
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of the central stellar mass and correspondingly to an underestima- 
tion of Mbh. 

In this work, we do not include a dark matter halo com- 
ponent in any of the mod els. Compared to the case of M87 
( iGebhardt & Thomas! I2009T) . the SMBH's sphere of influence in 
NGC 1332 is better-resolved by our SINFONI observations. The 
TVMbh degeneracy then becomes less severe as there are kine- 
matic constraints at several radii inside the sphere of influence. 
This kinematic data constrain the enclosed mass at each of those 
radii and help to disentangle the SMBH and the stellar contribution 
to the central potential. We decided that our measurements are suf- 
ficiently robust and it is not necessary to further investigate Mbh 
by including a dark halo component also for the following reasons. 
The kinematic data that we used in the modelling only extend out to 
r < 30 arcsec or equivalently r < 3 kpc. At this radius, we expect 
only a small fraction of the total mass to be in the dark matter halo. 
Empirically, the dark matter fraction found in early-type galaxies 
with a similar luminosity to NGC 1332 (Mb = — 20.5, as given 
by Hy perLeda) is at most 25% at a radius of 3 kpc dThomas et al.l 
2007). The mass decomposition for NGC 1332 in H09 suggests 
that the dark halo only takes ~10% of the total mass at 3 kpc. Fur- 
thermore, in our preferred model, the galaxy is decomposed into a 
bulge and a disc with each component having its own mass-to-light 
ratio. The bulge is well inside the region where the dark matter 
is expected to be unimportant. Its mass-to-light ratio is then 
not affected by the assumption about the halo. Neglecting the halo 
could bias T<j to be too high. However, the central stellar mass is 
dominated by the bulge so it does not suffer from any bias in 
and therefore the black hole mass is not affected by the exclusion 
of the dark halo. It is obvious from Run 2B and 2B* that Tt is 
independent of the variation in T^. 



7 SUMMARY AND DISCUSSION 

We presented SINFONI observations of the lenticular galaxy NGC 
1332 in the A"-band with the purpose to measure the mass of the 
SMBH. The sphere of influence of the SMBH (~0.76 arcsec) is 
resolved by the data (spatial resolution: FWHM ~ 0.14 arcsec). 
The kinematics, derived by fitting the CO bandheads, show a mod- 
erate rotation (v « 150 — 160 km s _1 at the outermost radius of 
SINFONI data) and a high velocity dispersion (~ 400 km s _1 in 
the centre). To complement our SINFONI data at larger radii, we 
utilised major-axis long-slit data of K96. Both datasets are con- 
sistent outside r > 2 arcsec; inside that radius the comparison is 
difficult due to the PSF difference. 

To determine the SMBH mass, we performed axisymmetric 
Schwarzschild modelling. The systematic uncertainties inherent in 
the dynamical models, which are due to assumptions on the inclina- 
tion, the density profile and the assumption of axisymmetry, were 
investigated. For the modelling input, two density models were 
constructed: a one-component model and a two-component model; 
the latter consists of a bulge and a disc. For each density model, 
we performed runs with and without the long-slit data. Using both 
SINFONI (r < 3.5 arcsec) and long-slit data (3.5 < r < 30 
arcsec) as the kinematic constraints, SMBH masses obtained using 
two different density models are consistent with each other. The 
inclusion of long-slit data appears to give significant constraints to 
the models and reduce the errors. The orientation of NGC 1332 is 
close to edge-on and the small uncertainties in the inclination as- 
sumption causes a negligible effect on the SMBH mass. The lack of 
obvious signs of triaxiality, e.g., isophotal twists or kinematics mis- 



alignments, justifies the assumption of axisymmetry. Mbh values 
measured in the four different quadrants are consistent with each 
other. 

Our preferred model is based on the two-component density 
profile, includes both SINFONI and long-slit data and adopts an 
inclination angle of 90°. We find a SMBH mass of Mbh=(1-45 ± 
0.20) x 10 9 M Q and a bulge mass-to-light ratio T 6 =7.08 ± 0.39 
in the R band. The disc mass-to-light ratio is not well-constrained 
by the data but falls within a range of 5.0 to 9.0 (also in i?-band). 
The comparison of our results to previous measurements and the 
implication of this A/bh for the A/bh-o" and the A/bh-£ relations 
are discussed in the following. 

7.1 Comparison with results from X-ray data 

Our results were preceeded by an X-ray analysis of NGC 1332 us- 
ing Chandra data (H09). They derive a most probable SMBH mass 
of 5.2±2:8 x 10 s M (5.4±|J x 10 8 M Q at the distance adopted 
in this work), and a mass-to-light ratio in J-band of 1.161q'j4 
(1.111q'i4 at our distance). As shown in Section [6] models with 
A/bh=5 x 10 8 A/q do not provide a proper fit to the kinematic 
data and they are separated from our preferred model by A\ 2 of at 
least 40. This is understandable since the SINFONI data give a high 
central velocity dispersion inside the sphere of influence and thus 
such a SMBH mass would be incommensurate. The A/bh value 
of H09 is obtained by using the A/bh-ct relation as the Bayesian 
prior. They are unable to rule out the possibility that the SMBH 
mass is an underestimate. The use of X-ray data alone, without the 
Mbh-o" prior, decreases their A/bh estimate even further although 
the upper limit is similar to our A/bh estimate (see their Fig. 8). 

Comparing our stellar mass profile with that of H09, we find 
that while the slope is very similar, our enclosed mass is twice 
as high at all radii where the stars dominate. The discrepancy 
is, however, not specific to H09; stellar mass profiles for NGC 
1332 derived from other X -ray analyses (Fuk azawa et al.l [20061 : 
iNagino & Matsushitd l2009h seem to be in agreement with H09. 
We find that lowering the inclination to 80° only has a negligible 
effect to the mass profile. From the orbital structure of our pre- 
ferred model we find that the rotation is only important at the outer 
radii. In Fig. [9] we plot v^/(u) along the major axis, where the ro- 
tation is highest. This quantity measures the importance of rotation 
with respect to random motion. At the inner radii, the significance 
of rotation is low and thus negligible for the mass estimation. The 
v$/{a) profile, however, rises with radius and peaks at a value of 
order unity at r ~ 20 arcsec. From this, we can expect that if the 
gas follows the stellar rotation, the X-ray-derived mass would be 
underestimated the most around this radius by roughly 25%. Rota- 
tion can therefore not be responsible for the systematic difference 
in the mass profile. The difference in mass is reflected by the dif- 
ference in the resulting mass-to-light ratios in our models and the 
models of H09. From the stellar population point of view, the low 
mass-to-light ratio measured by H09 can be obtained if this galaxy 
is young. However, the age of 4 Gyr as derived in H09 has a large 
upper error which is consistent with the old age found by Ba rr et all 
d2007l) (see subsectionlSTH. 

This is not the first ti me that a dynamically derived mass dif- 
fers f r om that of the X-ray dShen & Gebhardll2009l ; Ijohnson etail 
2009; iRomanowskv et alj 120091) . In those works, the invalidity 
of the assumption of hydrostatic equilibrium is suspected to 
contribute at least to some fraction of the discrepancy, as is 
also p ointed out bv lDiehl & Stalled j2007h and|ciotti & Pel legrinil 
(2004). Inflowing gas and the presence of non-thermal pressure 
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Figure 9. The significance of rotational velocity relative to the random mo- 
tion of the stars as a function of radius along the major axis, is the mean 
rotation in the azimuthal direction where + <r^ ) constitutes the second 
moment of the azimuthal velocity, (a) is the local mean velocity dispersion 
defined as (a)' 2 = (af, + ct| + c?)/3. The shade shows the area where all 
models within \u error fall and the line represents the average. 
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Figure 10. M Bii -cr (left) and M B h-L (rig ht) diagrams. NGC 133 2 is plot- 
ted as a square in each panel along with the Giiltekin et al. 1 2009) Mbh-v 
relation and Marconi&Hunt (2003) Mbh~L k relation. 



due to magnetic fields, microturbulence or cosmic rays are pos- 
sible reasons as to why the X-ray derived mass can be lower 
than the true mass jCiotti & Pellegrimll2004l : Ijohnson et al"1l2009l : 
IChurazov et alj|2008l) . In the case of NGC 1332, our dynamical 
modelling implies a much higher Mbh, which might indicate the 
failure of hydrostatic approximation, at least in the inner part. In the 
outer region, rotation can possibly account for a small part of the 
missing X-ray mass. All in all, the mass discrepancy probably in- 
volves a combination of systematic uncertainties in both methods, 
for which a detailed inspection is outside the scope of this paper. 
A situation similar to that o f NGC 1332 has been recently reported 
bv lShen & Gebhardtl d2009l) for NGC 4649. Their orbit-based mod- 
elling results in a larger Mbh and a ~ 70% higher mass profile than 
obtained from X-rays in H09. A larger sample of galaxies would in- 
deed be required to investigate whether this trend applies generally 
for these two methods. 



7.2 SMBH-bulge relation 

With our Mbh measurement, and our photometric decomposition, 
NGC 1332 is displaced from the standard relations in both Mbh-o" 
and Mbh-L d i agram s (Fig. [Tot . The galaxy is located above the 
iGultekin et ail ( 120091) M B n-a relation by 0.15-0.20 dex which is 
still within the intrinsic scatter. The velocity dispe rsion of 327.7 
km s - 1 was calculated using the definition of a e in Giilte kin" et al.l 
(2009). It was measured using the slit data within the effective ra- 
dius of the bulge of 8.39". Without the luminosity- weighting, the 
velocity dispersion drops to 319.2 km s _1 which is very close to 
the value found in HyperLeda. 

In the Mbh-L di agram, NGC 1332 is strikingly off of the 
iMarconi & Hunt! d2003l) relation. Our Mbh is one order of magni- 
tude higher than expected for its bulge luminosity. If the Mbh-L k 
(Lk is the bulge luminosity in the A'-band) were obeyed, it would 
result in a black hole mass of 1.37 x 10 s Mq which is highly ex- 
cluded in any of our runs. 

The a-L relation in the current black hole samples is known 
to be different from that in the SDSS sample, on which the distri- 
butions of L and a are based; the black hol e samples have larger 
a for a given L or smaller L for a given a dBernardi et al.ll2007l : 
iTundo et alj2007l:lLauer et al.ll2007l) . In this case NGC 1332 is not 
an excep tion. It is in fact a ra ther extreme outlier in the a-L di- 
agram of iBernardi elatl d2007l) for SMBH sample, i.e. the bias is 
stronger than expected from the SMBH sample. The Mbh-ct and 
the Mbh-L relations predict different SMBH masses and contra- 
dict each other by a factor of about seven. 

Provided that the bias in the a-L relation is just a selection 
effect, a question arises: which relation is the more fundamental 
one? Our result for NGC 1332 favours the A/bh-o" to be the more 
fundamental relation since the measured A1 B h for this galaxy turns 
out to fall much closer to the value predi cted by the Mbh-Q" r ela- 
tion. This is in line with the suggestion of lBernardi et al.l d2007l) . In 
addition, NGC 1332 is located slightly above the M B H-a relation, 
which makes it also consistent w ith the sugges tion that the relation 
curves upwards at the upper end (Wyithe 2006) or that the intrinsic 
scatter increases in this regime. 



ACKNOWLEDGEMENTS 

We thank the Paranal Observatory Team for support during the ob- 
servations. We are grateful to Koenraad Kuijken for sharing the 
long slit spectral data which we used for the dynamical mod- 
elling and to Ortwin Gerhard for useful discussions. We would 
also like to acknowledge the anonymous referee, whose comments 
have improved this paper. The research of P.E. is supported by the 
Deutsche Forschungsgemeinschaft through the Priority Programme 
1 177 'Galaxy Evolution' . Support for N.N. is provided by the Clus- 
ter of Excellence: 'Origin and Structure of the Universe'. 



REFERENCES 

Abuter R., Schreiber J., Eisenhauer F, Ott T, Horrobin M., Gille- 

sen S., 2006, New Astronomy Review, 50, 398 
Baes M., Dejonghe H., De Rijcke S., 2000, MNRAS, 318, 798 
Barr J. M., Bedregal A. G., Aragon-Salamanca A., Merrifield 

M. R., Bamford S. P., 2007, A&A, 470, 173 
Beifiori A., Sarzi M., Corsini E. M., Dalla Bonta E., Pizzella A., 

Coccato L., Bertola F, 2009, ApJ, 692, 856 
Bender R., 1990, A&A, 229, 441 



The central black hole mass ofNGC 1332 13 



Bender R., Moellenhoff C, 1987, A&A, 177, 71 
Bernardi M., et aL, 2006, AJ, 131, 2018 

Bernardi M., Sheth R. K., Tundo E., Hyde J. B., 2007, ApJ, 660, 
267 

Bertin E., Mellier Y., Radovich M., Missonnier G., Didelon P., 
Morin B., 2002, in D. A. Bohlender, D. Durand, & T. H. Hand- 
ley ed., Astronomical Data Analysis Software and Systems XI, 
Vol. 281 of ASP Conf. Ser., The TERAPIX Pipeline, p. 228 

Bonnet H, et aL, 2004, The Messenger, 1 17, 17 

Cappellari M., et aL, 2006, MNRAS, 366, 1 126 

Churazov E., Forman W., Vikhlinin A., Tremaine S., Gerhard O., 
Jones C, 2008, MNRAS, 388, 1062 

Ciotti L„ Pellegrini S., 2004, MNRAS, 350, 609 

de Souza R. E., Gadotti D. A., dos Anjos S., 2004, ApJS, 153, 41 1 

Di Matteo T, Springel V., Hernquist L., 2005, Nature, 433, 604 

Diehl S., Staffer T. S., 2007, ApJ, 668, 150 

Eisenhauer E, et al., 2003, in M. Iye & A. F. M. Moorwood ed., 
SPIE Conf. Ser. Vol. 4841 of Society of Photo-Optical Instru- 
mentation Engineers (SPIE) Conference Series, SINFONI - In- 
tegral field spectroscopy at 50 milli-arcsecond resolution with 
the ESO VLT. p. 1548 

Ferrarese L., Merritt D., 2000, ApJ, 539, L9 

Fukazawa Y., Botoya-Nonesa J. G, Pu J., Ohto A., Kawano N., 
2006, ApJ, 636, 698 

Gadotti D. A., 2008, MNRAS, 384, 420 

Gebhardt K„ et al., 2000a, ApJ, 539, L13 

Gebhardt K„ et al., 2000b, AJ, 119, 1157 

Gebhardt K„ et al., 2003, ApJ, 583, 92 

Gebhardt K., Thomas J., 2009, ApJ, 700, 1690 

Gerhard O. E„ 1993, MNRAS, 265, 213 

Granato G. L., De Zotti G., Silva L., Bressan A., Danese L., 2004, 
ApJ, 600, 580 

Grillo C, Gobat R., Lombardi M., Rosati P., 2009, A&A, 501, 
461 

Gultekin K., et al., 2009, ApJ, 698, 198 
Hiiring N., Rix H, 2004, ApJ, 604, L89 

Humphrey P. J., Buote D. A., Brighenti F, Gebhardt K., Mathews 

W. G., 2009, ApJ, 703, 1257 (H09) 
Johnson R., Chakrabarty D., O'Sullivan E., Raychaudhury S., 

2009, ApJ, 706, 980 
Kormendy J., Richstone D., 1995, ARAA, 33, 581 
Kuijken K., Fisher D., Merrifield M. R., 1996, MNRAS, 283, 543 

(K96) 

Kundu A., Whitmore B. C, 2001, AJ, 122, 1251 
Lauer T. R., et al., 2007, ApJ, 662, 808 
Magorrian J., 1999, MNRAS, 302, 530 
Magorrian J„ et al., 1998, AJ, 115, 2285 
Maraston C, 2005, MNRAS, 362, 799 
Marconi A., Hunt L. K„ 2003, ApJ, 589, L21 
Mei S., et al., 2005, ApJ, 625, 121 
Modigliani A., et al., 2007, arXiv:astro-ph/0701297 
Nagino R., Matsushita K., 2009, A&A, 501, 157 
Netzer H, 2003, ApJ, 583, L5 

Nowak N., Saglia R. P., Thomas J., Bender R., Davies R. I., Geb- 
hardt K„ 2008, MNRAS, 391, 1629 

Nowak N., Saglia R. P., Thomas L, Bender R., Pannella M., Geb- 
hardt K., Davies R. I., 2007, MNRAS, 379, 909 

Ogando R. L. C, Maia M. A. G, Pellegrini P. S., da Costa L. N., 
2008, AJ, 135, 2424 

Prugniel P., Heraudeau P., 1998, A&AS, 128, 299 

Romanowsky A. L, Strader L, Spitler L. R., Johnson R., Brodie 
J. P., Forbes D. A., Ponman T, 2009, AJ, 137, 4956 



Schreiber L, Thatte N, Eisenhauer F, Tecza M., Abuter R., Hor- 
robin M., 2004, in F. Ochsenbein, M. G. Allen, & D. Egret 
ed., Astronomical Data Analysis Software and Systems XIII, 
Vol. 314 of ASP Conf. Ser., Data Reduction Software for the 
VLT Integral Field Spectrometer SPIFFI. p. 380 

Schwarzschild M., 1979, ApJ, 232, 236 

Shen J., Gebhardt K., 2009, ApJ, submitted (arXiv:0910.4168) 
Shields G. A., Gebhardt K„ 2004, Bull. AAS, 36, 766 
Shields G. A., Salviander S., Bonning E. W., 2006, New Astron- 
omy Review, 50, 809 
Silge J. D., Gebhardt K., 2003, AJ, 125, 2809 
Silk L, Rees M. L, 1998, A&A, 331, LI 
Siopis C, et al., 2009, ApJ, 693, 946 

Somerville R. S., Hopkins P. F, Cox T. J., Robertson B. E., Hern- 
quist L., 2008, MNRAS, 391, 481 

Thomas L, Saglia R. P., Bender R., Thomas D., Gebhardt K., 
Magorrian L, Corsini E. M., Wegner G, 2007, MNRAS, 382, 
657 

Thomas L, Saglia R. P., Bender R., Thomas D., Gebhardt K., 
Magorrian L, Richstone D., 2004, MNRAS, 353, 391 

Tonry J. L., Dressier A., Blakeslee J. P., Ajhar E. A., Fletcher 
A. B., Luppino G. A., Metzger M. R., Moore C. B., 2001, ApJ, 
546, 681 

Trager S. C, Worthey G., Faber S. M., Burstein D., Gonzalez J. L, 

1998, ApJS, 116, 1 
Tremaine S., et al., 2002, ApJ, 574, 740 

Treu T, Auger M. W., Koopmans L. V. E., Gavazzi R., Marshall 
P. L, Bolton A. S., 2009, ApJ, submitted (arXiv:091 1.3392) 

Tundo E., Bernardi M., Hyde J. B., Sheth R. K., Pizzella A., 2007, 
ApJ, 663, 53 

van der Marel R. P., Franx M., 1993, ApJ, 407, 525 
Wyithe J. S. B., 2006, MNRAS, 365, 1082 
Wyithe J. S. B., Loeb A., 2003, ApJ, 595, 614 



APPENDIX A: SEEING CORRECTION IN THE 
DEPROJECTION 

Seeing due to atmospheric turbulence flattens the slope of the light 
profile of extended objects obtained by ground-based observations, 
especially at the innermost radii. When the surface brightness pro- 
file is deprojected without PSF correction, the resulting luminosity 
density will also be flatter than it actually is. Since light traces stars 
this will then lower the stellar mass contribution. The dynamical 
modelling only constrains the total mass, and so Mbu will be over- 
estimated to compensate for the decreased mass in stars. A priori, 
it is not known how much this effect alters the SMBH mass. 

We estimated the PSF effect using two one-component models 
described in Section |4~2l The first was the one-component model 
used in the main analysis, where we relied on the SINFONI pho- 
tometry (i.e., the collapsed datacube) for the innermost isophotes 
(hereafter referred to as the SINFONI dataset). The second one was 
the same model but we instead used the isophotes derived from the 
WFPC2 image to replace the SINFONI photometry at r < 0.5 arc- 
sec (HST dataset). We accounted for the seeing in the deprojection 
step described in Section |4~3l 

For the HST dataset, we generated a PSF for the Planetary 
Camera chip of WFPC2 with version 6.3 of the Tiny Tim software 
package0 using the location of the galaxy center and an appropri- 

2 http://www.stsci.edu/software/tinytim/tinytim.html 
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Figure Al. A double non-circular Gaussian fit to the HST PSF. The fit 
along the x and the y-axis of the detector are shown in the left and right 
panel respectively. The black dashed lines are the individual Gaussians and 
the red solid lines are the sum of the Gaussian components; diamonds are 
the actual PSF. 

ate K giant spectrum. This instrumental PSF of the WFPC2 image 
was then parametrised in the same way as that of the SINFONI im- 
age, i.e. as a non-circular double Gaussian function. Fig. lAll shows 
the fit to the HST PSF. The resulting FWHMs are 0.067 arcsec 
and 0.24 arcsec for the narrow and broad components respectively 
(FWHM X « FWFTM y for both components). 

The top panel of Fig. IA2I shows the luminosity density pro- 
files along the major and minor axes for both HST and SINFONI 
datasets. Both densities are almost identical at all radii along the 
major axis. There is only a small difference along the minor axis 
at 0.2 arcsec < r <0.8 arcsec, implying that the SINFONI density 
model is slightly more flattened at those radii. Due to the strong 
dust extinction in the WFPC2 image, the ellipticity in the inner part 
is not well constrained and is fixed to an approximate value. To- 
gether with the different shape of the HST and the SINFONI PSF, 
this prevents an exact match of SINFONI and HST density pro- 
files. However, the general agreement between the two profiles is 
very good considering that the inner profile was constructed using 
different datasets, which were deconvolved using different PSF im- 
ages. 

The bottom panel shows the reprojection of the luminosity 
density models of both datasets without seeing convolution, which 
should reflect the intrinsic surface brightness profiles. We see that 
both profiles agree well with each other in the inner part where the 
correction is most significant (r < 1 arcsec). In the outer part, the 
profiles overlap with each other, as is also the case in the luminos- 
ity density profiles. This is expected as both profiles use the same 
EMMI data at r > 4.5 arcsec and the effect of seeing is restricted 
only to the innermost isophotes. 

We performed the dynamical modelling using each of the lu- 
minosity density models from the two datasets, once with PSF cor- 
rection and once without. The modelling setups were identical to 
those in Section[5] Our parameter grid for each of the runs consisted 
of 20 trial values of Mbh ranging from 5 x 10 8 M o to 5 x 1O 9 M , 
each was paired with 20 different T values ranging from 1 to 10. 
To minimise computing time, we only used SINFONI kinematics. 
Mbh is sensitive to the change of the density profile and using only 
SINFONI data will put the least constraints to Mbh. Therefore in 
some sense these runs should show the largest possible change in 
Mbh due to the PSF inclusion. 

We list the best-fitting Mbh (average of four quadrants), 
marginalised over T, together with the averaged la errors in Ta- 
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Figure A2. The PSF-deconvolved luminosity density models (top panel). 
Thick and thin lines refer to density profiles along the major and minor 
axes respectively. The reprojection of these luminosity density models with- 
out seeing convolution resulted in the PSF-deconvolved surface brightness 
models (bottom panel). The SINFONI dataset (red dashed line) used SIN- 
FONI images out to 4.5 arcsec and the HST dataset (black solid line) used 
HST images to replace SINFONI data for the innermost isophotes (r < 0.5 
arcsec). For r > 4.5 arcsec, the NTT-EMMI image was used for both pro- 
files. 



Table Al. SMBH masses obtained for different modelling runs with differ- 
ent luminosity models. The given errors are the average of \a errors from 
the four quadrants. All A^bh are given in units of 10 9 AfQ. 



Luminosity model A/bh 



HST without PSF correction 


1.56 


± 


0.24 


SINFONI without PSF correction 


1.92 


± 


0.24 


HST with PSF correction 


1.27 


± 


0.22 


SINFONI with PSF correction 


1.15 


± 


0.24 



ble IA1I For comparison, we rewrite the results of run 1A in Table 
[T] as the model "SINFONI with PSF correction". When the see- 
ing effect on the photometry is taken into account, A/bh decreases. 
This is expected, as the PSF deconvolution steepens the slope of the 
surface brightness and the luminosity density profile, giving more 
mass to the stars. The mass-to-light ratio T increases only slightly, 
from ~ 8.6 for models without PSF correction to ~ 8.9 when ei- 
ther of the PSFs is included. The photometric difference (between 
the HST and the SINFONI image) in the innermost region does not 
seem to affect T. 

For the runs where the PSF is not included, the SINFONI 
Mbh is higher than that of HST. When the PSF is included, the 
SINFONI A/bh decreases more dramatically (see Table lAlt . These 
changes are expected, as the HST PSF is narrower than the SIN- 
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FONI PSF. When the HST data are used, the change of M BH due 
to the PSF is still within the lcr error, while for SINFONI data the 
change is larger, but still less than a factor of two. Although this ef- 
fect would not apply equally to all galaxies and observation modes, 
our exercise seems to suggest that the PSF can be considered a 
relatively minor issue for HST photometry. However, for ground- 
based observations which in general have broader PSFs, more care 
is needed. 

Comparing Af B H obtained using the HST and the SINFONI 
datasets after PSF deconvolution, we see that both masses are con- 
sistent and lie well within their la errors. This result is reassuring 
considering the uncertainties in the SINFONI PSF. 



